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(54) Method of and apparatus for digital radio signal reception 



(57) A receiver is described, for instance for GSM 
systems, wherein a plurality of digital signals each com- 
prising a training sequence (Xp (t)) usable to generate 
an estimation of the transmission channel are received 
in diversity (A.1 ... A.N) in the form of a plurality of repli- 
cas (10.1 ... 10. N) each comprising a respective replica 
of the training sequence. The received digital signal is 
subjected to a delaying action thereby obtaining a plu- 
rality of versions of such signal, each comprising a 
respective set of signal replicas. Each version is sub- 
jected to a respective filtering action (23.1 ... 23.N) inde- 
pendently from the other versions, by employing a 
respective first set of filtering coefficients (w*^ ... w* N1 
... w* N j) obtained starting from a respective initial set of 
filtering coefficients. This respective initial set is 
obtained by performing, for each replica of the received 
signal, a filtering action exploiting a respective second 
set of filtering coefficients identified as a signal (u (t)) 
able to generate, by convolution (u (t) * Xp (t)) with the 
training sequence, a unitary function 8 (t) on a given 
time window, preferably on one bit. The filtering coeffi- 
cients of the first set are derived starting from said initial 
set of filtering coefficients through an MMSE adapta- 
tion, preferably with a initial rapid convergence phase of 
the RLS type followed by an optimisation phase of the 
LMS type. 
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Description 

The present invention relates to digital communication and more specifically it concerns a method and an appara- 
tus for the reception of digital radio signals, according to the preamble of claim 1 . 

5 The invention has been developed with particular attention to its possible use within mobile radio communication 
systems such as the systems known as GSM and IS-54. In any case, the invention can be applied in any contest 
wherein the functional reception diagram applied is substantially similar to or identical with the diagram shown in Fig. 1 . 

To illustrate how the invention can be utilised in an existing system, the diagram of a conventional, single antenna, 
GSM receiver is presented in Rg. 1. In this diagram reference 10 indicates a line whereon a received baseband digital 

w signal is present. In addition to an actual useful signal component, the signal comprises a training signal component 
such as a preamble or a so-called "midamble", essentially comprising a string of binary symbols Xp (t) which is assumed 
to be known in compliance with the GSM standard. The signal present on line 10 is obtained from the radio signals 
received by a antenna through a processing in a module indicated by reference 20: there, the antenna signal is, for 
instance, converted into baseband and subjected to the conditioning operations currently used in the art in view of Hs 

15 demodulation. 

The signal on line 1 0 is split by a demultiplexing block (not shown in the drawing) between two branches 1 1 and 1 2, 
intended to convey the training signal or the information signal, respectively. It should be stated that, as is well known 
to the person skilled in the art, such splitting does not necessarily correspond to an actual routing over two different 
physical channels, since it can be effected in a virtual manner by means of different processing operations. 

20 The processing performed by the first branch 1 1 is aimed at obtaining an estimate of the channel impulse response 
whereon the transmitted signal has propagated. This estimated response can be obtained by analysing how the training 
signal is affected by the channel. The aforesaid estimate is usually performed by the correlation (or the matched filter- 
ing) effected in a module indicated by 14. In block 15, the convolution is calculated between the windowed estimate of 
the channel impulse response and a set of possible transmitted baseband signals S' p (t) (over one bit period) to obtain 

25 the signal estimates x (t). These estimates are fed to a processing module 16 where the signal routed over the branch 
12 arrives after a possible filtering effected in filter 17. This filter has a impulse response equal to the windowed ambi- 
guity function (the ambiguity function being, as known, the auto-correlation of the training sequence), i. e., a response 
given by [x p (t) • x p (T - 1)] h w (t) , where rv, (t) is a window function. This method has been described by R. Steel, 
"Mobile Radio Communications", New York 1992, Chapter 6. In processing module 16 a measurement of the "distance" 

so (incremental metrics for one bit interval) between the generated sequences and the actual received data is performed. 
The incremental metrics calculated in block 16 are fed through line 1 6 to a Viterbi processor (known in the art), included 
in block 19, where the new metrics for each state are established, as occurs in GSM transmission systems. The Viterbi 
processor is followed in cascade arrangement by a differential decoder (equally known in the art) which emits the output 
data stream. Essentially, the received signal on line 10 is submitted to a processing operation that can ideally be seen 

35 as a complementary and opposite action to the one effected by the transmission channel. 

From the analysis of the reference diagram in Fig. 1 , it is clear that the more articulated is the set of alterations the 
transmitted signal - and also the training sequence - may have undergone during transmission over the cannel, the 
more complex and onerous is the set of processing operations performed in elements 14, 15, 16 and 19. 

In particular, in mobile radio systems (at least for base stations, but the use of this technique is being extended also 

40 to mobile terminals) the use of diversity reception techniques based on the use of a plurality of N receiving antenna has 
become widespread. The signal received by an array of antenna of this kind in reality comprises N replicas of the same 
starting signal, which replicas are received by the aforesaid N antennas in a different way (for example due to a different 
distribution of the echoes, etc.). 

The invention exploits this multiplicity of antennas to develop a more robust receiver structure which enhances the 

45 communication link quality. The processing of the individual symbol at the receiving side entails the analysis of a certain 
number (for example, M) of symbols received successively. A receiver architecture with N antennas should, for the 
same propagation environment, consider the analysis of NxM symbols. Recovering the transmitted signal in case of the 
multi-input system essentially entails inverting a system matrix N x j, where N is the number of antenna and j the 
number of time instants 0- a- the number of consecutive snapshots) considered necessary to faithfully reproduce the 

so transmitted signal. Direct inversion of this matrix can, if done without due attention, lead to noise amplification and insta- 
bility. Moreover, it can be rather onerous in terms of time and hardware required, and it hardly appears practical for real- 
time processing of the received signals, as is required in case, for example, of voice signals. 

The object of the present invention therefore is to provide a solution that, though similar to the general diagram 
shown in Fig. 1 , does not give rise to the drawbacks described above, further allowing a greater resolution in the per- 

55 formance of the auto-correlation function of the training signal Xp (t). 

According to the present invention, this object is attained thanks to a method and a system having the characteris- 
tics set out specifically in the claims that follow. 

The invention is primarily based on the fact that in a field of use such a the GSM mobile radio system, an echo over 
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the individual channels of diversity reception (N receiving antenna) is essentially discretised by sampling the received 
signal, and therefore the matrix that identifies the behaviour of the channel (i. e. the matrix to be estimated and inverted 
at the receiving side) actually comprises a number of elements with null or nearly null values. 

The invention shall now be described, purely by way of non limiting example, with reference to the accompanying 
5 drawings, wherein: 

- Fig. 1 schematically shows a receiver architecture for GSM systems, that has already been extensively described 
above; 

Fig. 2 shows, in block diagram form, the structure of a receiver operating according to the invention; 
10 - Fig. 3 is a diagram illustrating the performance of the system according to the invention in terms of auto-correlation 
of the training signal; and 

- Fig. 4 is a diagram illustrating the performance of the invention in a given transmission channel; 

- Fig. 5 illustrates a variant of the invention. 

15 The diagram in Fig. 2 illustrates an extended general layout block diagram which is similar to the one shown in Rg. 
1 and makes use of a plurality of input antenna A.1 ... A.N, connected according to a typical diversity reception scheme. 
Like elements in Figs. 1 and 2 are indicated by like references. 

The single input line 10 shown in Fig. 1 corresponds in Rg. 2 to N input lines 10.1 ... 10.N, each conveying a 
respective replica of the signal received, which can be expressed, in general, in the form of a complex signal, thus a 

20 signal with a real part and an imaginary (or quadrature) part Each of these replicas is derived from a respective input 
antenna A.1 ... A.N through a conversion from radio frequency to intermediate frequency followed by a conversion from 
the intermediate frequency to the base band. Such conversions are individually performed in a known manner on the 
N signal streams coming from the N antenna in blocks 20a and 20b corresponding, as a set, to block 20 in Figure 1 . 
Reference 26 indicates a demultiplexing unit that splits the information signal and the known midamble into different 

25 streams in the same manner as mentioned in connection with Fig. 1 , the only difference being that the operation per- 
formed by demultiplexer 26 correspond to the splitting of N parallel paths. 

Blocks 16 and 19 (the latter being split in the Figure into Vrterbi decoder 19a and differential decoder 19b) perform 
the same operations as the homonymous blocks in Fig. 1 . 

The whole of the circuits enclosed in dashed line block 22 perform the "channel estimation" of the sampled base- 

30 band received signals 10.1 ... 10.N according to the invention. The estimated channel impulse responses from the out- 
puts of FIR (Finite Impulse Response) fitters 27.1 ... 27.N are achieved by computing the convolution between the input 
training signals Xp (t) (provided by block 29) and the impulse responses u (t) of the Fl R filters. These filters can be imple- 
mented as transversal filters. The impulse response u (t) is determined by using the assumption that the convolution Xp 
(t) * u (t) is a quasi- Dirac impulse. By sampling impulse response u (t) in block 28, the coefficients of the FIR filters are 

as obtained. The selection of the appropriate impulse responses u (t) depends of the propagation channel. The outputs of 
filters 27.1 ... 27.N are processed in a block 25 where the initial weights for the adaptation algorithm performed in block 
24 are selected. The adaptation algorithm is advantageously an MMSE (Minimum Mean Square Error) algorithm and 
can be implemented as an RLS (Recursive Least Square) or an LMS (Least Mean Square) algorithm. 

Block 23 (beam-forming network) comprises a bank of j filters individually shown at 23.1 ... 23.j. Each filter operates 

40 on a corresponding version of the input signal, a first version (processed by fitter 23.1) corresponding to the signal 
received on the lines 10.1 ... 10.N whilst the other versions correspond to the same signal successively delayed by 
delay intervals T, 2T ... (j - 1) T as an effect of the passage through corresponding delay lines 230. Interval T can cor- 
respond, for instance, to the bit period used for transmission (conventional equally spaced equaliser). However, frac- 
tionally spaced or non-equally spaced equalisers can also be used. 

45 More particularly, in filter 23.1 weighting coefficients w*^ , w* 21 .... w* N1 are applied to the replicas of the received 
signal. A similar processing occurs in filter 23.2, with the application, to the various signal replicas delayed by an interval 
T, of coefficients w* t 2 , w* 13 ... w* N2 , and so on to filter 23,j, where coefficients w*!], w* 2 j ... w* Nj are applied to the vari- 
ous signal replicas delayed by an interval Q - 1) T. 

The aforesaid weighting coefficients, which can be ideally represented as a matrix with N rows and j columns, are 

so computed by block 24 that applies a mechanism of estimation with adaptation which will be better described further on. 
Regardless of the possible presentation of the coefficients w* u ... w* N j in matrix form, it will be appreciated that the 
set of filters 23.1 ... 23.j effect, on the various versions of the signal present on lines 10.1 ... 1 0.N, a processing action 
that does not have a matrix nature in each of the j time stages: indeed, each of the various versions of the signal suc- 
cessively delayed by an interval T is subjected to an independent filtering operation (see e. g. Fuhal and Bonek, "Space- 

55 time decomposition: Exploiting the full information of a training sequence for an adaptive array", Electronics Letters, vol . 
32, no. 21, pp. 1938 - 1939, October 1996). 

This choice has proved advantageous with respect to the general solution such a that mentioned above, where the 
channel estimation and equalization or beam forming of the received signal entail a matrix estimation and inversion 
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operation, with the drawbacks analysed above. Differently from that known solution, the embodiment shown in Fig. 2 
provides for a channel estimation operation performed independently on each version of the received signal, consid- 
ered in the set of the replicas it is composed of: this allows taking into account the effects of the channel (typically due 
to echo phenomena) on that version. This results, on one hand, into a simplification of the processing operations and, 
5 on the other hand, into the possibility of avoiding negative phenomena (error propagation, noise amplification) linked to 
the specific nature of the received signal (for instance a GSM signal, in the exemplary embodiment considered). 

The following mathematical description of the method deals with a very important aspect of the invention, that is 
the improved channel estimation in the initial step of the weight adaptation. In particular, the coefficients of the impulse 
response u (t) in FIR filters 27.1 to 27.N are chosen in such a way as to meet the relation 

10 

x p (t)'u(t) = 8(t) (D 

where 5 is the Dirac delta. 

Impulse response u (t) is provided to block 30, which actually effects the convolution defined in equation (1). The 
is output of block 30 is multiplied in multiplier 32 by a window function (t) to produce a windowed ambiguity function ^ 
(t)- 

The output signal of beam forming network 23 is fed to a processing unit 31 which provides on a first output 31 a 
the received training signal Xp (t) affected by a scalar channel response a (t), i. e., 

x p (t) = x p (t)*a(t) (2) 

By the convolution of that training signal by the impulse response u (t) in block 38 and multiplication by the same window 
function h* (t) in a multiplier 34 an estimate a* (t) for a given delay is obtained, according to relation: 

25 [x p (t)*a(t)*u(t)]h w (t) (3) 

As known to those skilled in the art. the convolution performed in block 30 is the more delta-like the more accurate 
is the estimate of the channel response. However, this desired property must be compromised to a same extent in 
actual real-time communication systems. 
30 A second oulput 31b of processing unit 31 conveys the information signal x (t, x) affected by the same impulse 
channel response a (t). The convolution between such signal and the windowed ambiguity function (t) is computed 
in block 1 7, that gives the following output: 



x(tt)*a(t)*6 w (t) (4) 

35 

The output of block 1 7 is used in block 16 for calculating the incremental metrics. 

Block 15 performs the convolution between the windowed channel response aw (t) and the baseband signal gen- 
erated in modulator 1 5'. In this way all possible signals are created. They are used in the block 1 6 for obtaining all incre- 
mental metrics as in Fig. 1 . 

40 The determination of the weights in block 24 is based on the minimisation of the mean square error between the 
received training signal on 31a and the original training signal Xp (t) from block 29, already known to the receiver. The 
subtraction between the two signals is performed in 42. 

The receiver structure shown in Fig. 2 is adapted to be used with particular advantage, for instance, within a GSM 
transmission scheme wherein a midamble of 26 bits is surrounded by 58 information bits. 

45 In this case, the midamble sequence a n with n=1 ... 26 (physical values) is characterised by a particular auto-cor- 
relation property which allows channel estimation with a 5 bit delay, i. e. an estimation based on relation: 



IfiZ^k+n^n n = " 4 0....4 (5) 

50 k=5 



However, whilst this property expresses a Dirac delta on one bit, when the modulation is applied, the auto-correla- 
tion of the modulated sequence Xp (t) expressed as 

55 

RM=fXp(t)*X p (T-t)dt (6) 

is smeared out over 4 bits, with a consequent loss of resolution. 



4 



10 



EP0 892 504 A2 

* In the situation just described it is possible and particularly advantageous to use a method of filtering the training 
sequence that originates an auto-correlation function having a delta-like behaviour on a smaller interval, in particular, 
but not exclusively, of one bit. 

The calculation of the improved initial channel estimate (IICE) response is obtained by using the matrix represen- 
tation of the sampled signals Xp (n) and u (n) . For instance the output of a matched filter like block 1 4 in Fig. 1 is given by 



Rmf(")= Z x p (n-k).x p H (k) (7) 
kol 

where Xp (n) is the known transmitted training sequence, length (x) indicate the sequence length and the superscript H 
indicates hermitian conjugate. The response of the IICE filter, u (n), on the signal Xp (n) is given by 

15 length (x) 

Rmf(")= Z x p (n-k).x p H (k) (*) 

k-1 



20 The auto-correlation function determined by relation (7) has a time duration equal to twice the time duration of the 
GMSK impulse response (i. e., 4 bit times). We can perform a better auto-correlation function by using an impulse 
response u (n) which satisfies the following matrix equation 



25 



=8 



(9) 



where 



30 



V 



*16 



*11 



(10) 



35 



The equation system expressed above needs not necessarily be square, since the number of constraints can be 
greater, less than, or equal to the number of variables in u (t). 

The inversion of this system is performed by using singular value decomposition and the related solution can be 
40 expressed in the form 



u(t) = X" 1 -5 



(11) 



45 The auto-correlation of this filter and that of the conventional GSM matched filter are comparatively shown in Fig. 
4 where the number of bits is plotted on the x-axis and the scale of the y-axis is a scale normalised so as to attribute a 
unit value to the auto-correlation peak. The solid line curve represents the conventional estimation, whereas the dashed 
line represents the behaviour of the innovative filter described above. It will be noted in particular that the latter origi- 
nates an auto-correlation on a single bit with virtually total suppression of the auto-correlation value for a band of a few 

so bits respectively before and after the central auto-correlation peak. 

Channel estimation using the conventional approach and the approach proposed herein can be expressed respec- 
tively by the two equations that follow: 



hest_oid(*)= Jy»xp VQdt 02) 
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hesUiewM* Jy»«M* ( 13 ) 



5 

where y (t) = x p (t) * h (t) , h (t) representing the pulse response of the transmission channel, and x* p is the conjugate 
complex of Xp(t). 

Since the correlation between u (t) and Xp (t) approximates the behaviour of an ideal Dirac delta function better than 
equation (12) reported above, by using equation (13) rather than equation (12) a better channel estimation resolution 
10 is obtained. 

The coefficients of u (t) can be further multiplied, if deemed appropriate, by fixed coefficients to form further spatial 
filtering such as a Chebychev filtering. 

Note that a filter like filter 17 in Fig. 1 is now no longer necessary (in the exemplary embodiment of application to 
GSM) as the matched filter is not used in the channel estimation. Such filter may however be present in case of appli- 
75 cation to other systems. 

From the above it follows that the solution according to the invention allows obtaining a receiver wherein channel 
interference and the non-linear effects of the modulation, such a GMSK modulation, are treated separately. The 
receiver is capable of exploiting the spatial structure present both in the signal and in the interference. On the basis of 
these premises it is possible to obtain a receiver structure operating both in angle (diversity reception) and in time in 

20 order to contrast both intersymbol interference (ISI) and interference from different transmitters. 

In particular, the diagram in Fig. 2 shows the use of an MMSE (Minimum Mean Square Error) equaliser to mitigate 
the channel defects separately, followed by a Vrterbi decoder for the equalization of the GMSK signal (demodulation). 
The signal coming from the MMSE equaliser is equalised and phase aligned. The receiver can be utilised so a to incor- 
porate different diversity reception schemes, with a solution that, at least at present, is particularly advantageous for use 

25 in the base stations of mobile radio systems. The same structure can in any case also be applied to mobile receivers 
equipped with multiple antennas for diversity reception. 

The salient characteristics of the solution shown - purely by way of non limiting example - with reference to the cur- 
rently preferred embodiment of the invention can be briefly described in the following terms. 

The signal received by the N antenna elements A.1 A.N is converted into base band, and the digitised signal 

30 thus obtained is demultiplexed (block 26) thereby separating the training sequence (typically a midamble) from the data 
sequence. The filtering technique described above is applied to the N sequences thus obtained thereby attaining the 
identification of the relevant echoes. The values corresponding to the echoes identified (in practice the complex conju- 
gates thereof) and their range are used as initial values of the MMSE algorithm. The subsequent iterations of the related 
scheme (typically RLS) are performed on each identified echo using the training sequence. In particular this adaptation 

35 can be applied to the first and the last bit of the training sequence thereby to avoid miscalculating interfering signals. 
These outputs are further adapted for each echo by using an LMS algorithm which, if necessary, can be used to invert 
the entire system (all echoes). 

This solution further allows better echo identification and a reduced white noise gain since only the echoes are 
processed. A good initial estimation of the weights is also available, based only on the signal, together with a rapid con- 

40 vergence of the RLS algorithm, since the latter operates only on N weights for each echo, with the consequent possi- 
bility of performing the related processing operations in parallel on the various channels: thus, both a reduction of the 
calculation burden and a quicker calculation are obtained just thanks to the possibility of operating in parallel. It is pos- 
sible to think of adopting an LMS type iteration in order to obtain an overall optimal scheme starting from the sub-opti- 
mal echo-by-echo approach. As stated above, the possibility of using bits at the beginning and at the end of the training 

45 sequence leads to a very significant reduction in the possibility of an incorrect detection of interfering signals coming 
from other transmitters. 

The advantages of the above described invention can be summarised as follows: 

1) The initial weight adaptation helps coming closer to the global minimum of the MMSE function, and therefore 
so improves the convergence speed of the algorithm used; 

2) Block 22 actually provides for a suppression of unwanted signals from the other users (co-channel interference 
suppression): 

3) Using the space-time equaliser results in adjusting the signal from the antenna elements so as to give maximum 
signal-to-noise ratio by using a maximal-ratio combining strategy, since the adaptation of different stages is carried 

55 out in parallel, 

4) The improved initial channel estimate suggests the possibility of using a maximum likelihood sequence estima- 
tion as the optimal strategy for dealing with intersymbol interference, since the co-channel interference has been 
suppressed in block 22: In such conditions the block diagram in Fig. 2 may be adapted to the so-called blind equal- 
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isation and simplified as shown in Fig. 5, where an MMSE equaliser 100 receiving the data signal through line 31b 
and the weights present at the output of block 24 replaces the whole of blocks 17, 32, 33, 34, 15', 15, 16 in Fig. 2 

The preferred embodiment of the invention illustrated in Fig. 2 has the additional advantage of an improved syn- 

5 chronisation property. The synchronisation strategy is based on the oversampling of the received signal. The received 
signals are sampled with a sample rate P times the bit rate. Thus P different data sequences are obtained. In this way 
the resolution of the improved initial auto-correlation function is better and the clock synchronisation is improved. There- 
fore a coarse synchronisation can be obtained by sampling the received signal at the bit rate regarding to the embodi- 
ment in Fig. 2. A fine synchronisation will be achieved by using the cross-correlation property of the eight training 

10 sequences defined in the GSM standard, but in this case the received signals should be oversampled. A proper thresh- 
old level helps selecting the auto-correlation function of the desired user and regarding to this a fine sysnchronisation 
will be established. Essentially, this information is obtained from block 25 in Fig. 2. 

The improvement in system performance obtained by the use of the invention can be appreciated by considering 
Fig. 4, which shows the bit error rate BER versus ratio Eb/No (Eb = energy per bit, No = spectral noise density) for a 

75 GSM channel when using the invention (dashed line) and with a conventional receiver. 

Even if the invention ha been described with particular reference to the use in connection with GSM system, it can 
be applied without undue experimentation also to other systems in which the transmitted signals include a training 
sequence known a-priori. In particular, the adaptation module 18 is capable of generating the set of filtering coefficients 
also in multi-user systems, such as SDMA (Space Division Multiple Access) systems, in which two or more users within 

20 the same ceil are allowed to use the same physical channel (time - frequency). In this case, it is necessary that said 
users use training sequences Xp with sufficient orthogonal characteristics allowing the algorithm to distinguish between 
echoes originated by the different users. 

Of course, without altering the principle of the invention, the implementing details may vary widely with respect to 
what has been described and illustrated, without thereby departing from the scope of the present invention. 

25 

Claims 

1. Method of receiving digital signals comprising a training sequence Xp (t) usable .to generate an estimate of the 
transmission channel, wherein the digital signal is received in diversity as a plurality of signal replicas (10.1 ... 10.N) 

30 each comprising a respective replica Xp (t) * a (t) of said training sequence, the method being characterised in that 
it comprises the following operations: 

- generating a plurality of versions of the received digital signal separated by a given delay interval (T) , each ver- 
sion comprising a respective set of signal replicas; 

35 - subjecting each of said versions of the received digital signal to a respective filtering action (23.1 ... 23.j) inde- 
pendently of the other versions; said respective filtering action being performed on each of said versions by 
applying to the respective set of signal replicas a respective first set of filtering coefficients (w*^ ...w* N1; ... ; 
w*fj ... w* N p obtained starting from a respective initial set of filtering coefficients; 

- obtaining said respective initial set of filtering coefficients by subjecting the respective version of the received 
40 digital signal to a respective second filtering action (27.1 ... 27.j), performed independently of the respective 

second filtering actions effected on the other versions of the received digital signal; each of said respective sec- 
ond filtering actions being performed with a respective second set of filtering coefficients identified as a signal 
(u (t)) able to generate, by convolution (u (t) * Xp (t)) with said training sequence, a unitary function on a given 
time slot. 

45 

2. Method as claimed in claim 1 , characterised in that said respective second set of filtering coefficients (u (t)) are 
identified a a signal able to generate by convolution (u (t)) * x p (t)) with said training sequence a unitary function 
on one bit. 

so 3. Method as claimed in claim 1 or 2, characterised in that said respective first set of filtering coefficients is generated 
starting from said respective initial set of filtering coefficients by means of a Minimum Mean Square Error (MMSE) 
adaptation. 

4. Method as claimed in claim 3, characterised in that said adaptation is chosen in the group constituted by a Recur- 
55 sive Least Square (RLS) adaptation and by a Least Mean Square (LMS) adaptation. 

5. Method as claimed in claim 4, characterised in that said minimum mean square error adaptation comprises an ini- 
tial rapid convergence phase of the recursive least square (RLS) type followed by a phase of the least minimum 
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square (LMS) type. 



6. Method of receiving digital radio signals including a training sequence Xp (t) usable to generate an estimate of the 
transmission channel, wherein the digital signal is received in diversity as a plurality of signal replicas (10.1 ... 10.N) 

5 each comprising a respective replica Xp (t) * a (t) of said training sequence, the method being characterised in that 
it comprises the operation of performing an initial estimate of the transmission channel by subjecting each replica 
of the training sequence to a respective filtering action (27.1 ... 27.N), performed independently of filtering actions 
effected on the other replicas of the training sequence; each of said filtering actions being performed with a respec- 
tive set of filtering coefficients identified as a signal (u (t)) able to generate, by convolution (u (t) * Xp (t)) with said 

10 training sequence, a unitary function on a given time slot. 

7. Apparatus for the reception of digital signals comprising a training sequence (Xp (t)) usable to generate an estimate 
of the transmission channel, the apparatus comprising a plurality of antenna operating in diversity, so that the digital 
signal is received in the form of a plurality of signal replicas (10.1 , ... 10.N) comprising each a respective replica (Xp 

is (t) * a (t)) of said training sequence, the apparatus being characterised in that it further comprises 

- delay means for generating a plurality of versions of the received digital signal separated by a given delay inter- 
val (T), each version comprising a respective set of signal replicas, 

- first filtering means (23.1 , 23.2 ... 23 J) for subjecting each of said versions of the digital signal received to a 
20 respective filtering action independently from the other versions; said respective second f iKering action being 

performed on each of said versions by applying to the respective set of signal replicas a respective second set 
of filtering coefficients (w*^ ... w* N1 ...w* 11f w*^ ... w* Nj ) obtained starting from a respective initial set of filtering 
coefficients determined, 

- second filtering means (27.1, 27,2 .... 27.N) for obtaining said respective initial set of filtering coefficients by 
25 subjecting the respective replica of said training sequence to a respective second filtering action performed 

independently from the respective second filtering actions performed on the other respective replica of said 
training sequence; each of these respective second filtering actions being performed with a respective second 
set of filtering coefficient identified as a signal (u (t)) able to generate by convolution (u (t) * Xp (t)) with said 
training sequence a unitary function on a given time slot. 
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8. Apparatus as claimed in claim 7, characterised in that said respective set (second set) of filtering coefficients (u (t)) 
is identified as a signal able to generate by convolution (u (t) * Xp (t)) with said training sequence a unitary function 
over a preset interval. 

35 9. Apparatus as claimed in claim 8, characterised in that it comprises an adaptation module (24) which calculates said 
respective first set of filtering coefficients starting from said respective initial set of filtering coefficients by means of 
an adaptation of the minimum mean square error type. 

10. Apparatus as claimed in claim 9, characterised in that said adaptation module (24) operates according to an algo- 
40 rithm chosen in the group comprising recursive least square adaptation and least mean square adaptation. 

11. Apparatus as claimed in claim 10, characterised in that said adaptation module (24) operates with an initial rapid 
convergence phase of the recursive least square type followed by a phase of the least mean square type. 

45 1 2. Apparatus as claimed in any of claims 9 to 1 1 , characterised in that said adaptation module (24) is arranged to gen- 
erate a respective first set of filtering coefficients for each of a plurality of users within of a same cell of a digital radio 
communication system in which said users use the same physical channel (time - frequency channel). 

13. Apparatus for the reception of digital signals comprising a training sequence (Xp (t)) usable to generate an estimate 
so of the transmission channel, the apparatus comprising a plurality of antenna operating in diversity, so that the digital 
signal is received in the form of a plurality of signal replicas (10.1 , ... 1 0.N) comprising each a respective replica (Xp 
(t) * a (t)) of said training sequence, the apparatus being characterised in that it further comprises filtering means 
(27.1 , 27,2 ... 27.N) for performing an initial estimate of the transmission channel by subjecting a respective replica 
of said training sequence to a filtering action performed independently from the filtering actions performed on the 
55 other replicas of said training sequence; each of these filtering actions being performed with a respective set of fil- 
tering coefficients identified a a signal (u (t)) able to generate by convolution (u (t) * Xp (t)) with said training 
sequence a unitary function on a given time slot. 
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